Thermal Radiation Prom Carbon Nanotube in Terahertz Range 
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The thermal radiation from an isolated finite-length carbon nanotube (CNT) is theoretically 
investigated both in near- and far-field zones. The formation of the discrete spectrum in metallic 
CNTs in the terahertz range is demonstrated due to the refiection of strongly slowed-down surface- 
plasmon modes from CNT ends. The effect does not appear in semiconductor CNTs. The concept 
of CNT as a thermal nanoantenna is proposed. 

PACS numbers; 78.67.Ch, 44.40.+a 



Carbon nanotubes — quasi-one-dimensional carbon 
molecules — have found numerous real and potential ap- 
plications as building blocks of nanoelectronic circuits [Ij 
and nanoscale optical elements Among others, the 
idea of CNT-based optical devices enabling control and 
enhancement of radiation efficiency on the nanoscale, i.e. 
nanoscale antennas for infrared and visible light, is ac- 
tively discussed [1, 0, H, @] . Noise properties and opera- 
tional limits of such devices are substantially determined 
by the thermal fluctuations of electromagnetic field. Fun- 
damental interest to the thermal radiation is also dictated 
by the ability of nanostructures to change the photonic 
local density of states CLDOS) , i.e. the electromagnetic 
vacuum energy 0, 0, S]- The effect has been observed 
in microcavities, photonic crystals, nanoparticles in the 
vicinity of surface-plasmon resonances, etc. 0]. Thus, as 
the electromagnetic fluctuations are defined by the pho- 
tonic LDOS, investigation of the thermal radiation is ex- 
pected to bring new opportunities for the reconstruction 
of photonic LDOS in the presence of nanostructures. In 
turn, the photonic LDOS is a key physical factor defining 
a set of well-known quantum-electrodynamical ejfects: 
the Purcell effect fs'], the Casimir-Lifshitz forces 9], the 
electromagnetic friction [2], etc. All aforesaid stimulates 
investigation of thermal radiation on nanoscale. 

During last decade, there has been considerable inter- 
est to optical properties of surface-plasmon structures 
[13, El- In 

particular, spherical gold particles in the 
vicinity of the plasmon resonance have been shown to be 
an effective nanoantennas [13] • Recently, the strong co- 
herent coupling between individual optical emitters and 
guided plasmon excitations has been predicted [isj intro- 
ducing thus quantum optics of nanoplasmonic structures. 
Thermal radiation in systems with surface plasmons is in- 
fluenced by the near-field effects and is known to be con- 
siderably different from the black-body radiation 0,113] ■ 
Earlier theoretical studies of CNTs showed the existence 
of low-frequency plasmon branches [Til . and the for- 
mation in CNTs of strongly slowed-down electromagnetic 
surface waves jl8l] . Such waves define pronounced Pur- 
cell effect in CNTs [13] and potentiality of CNTs as 
Cherenkov-type emitters [i^l- Geometrical resonances 
— standing surface waves excited due to the strong re- 



flection from the tips of finite-length CNTs — qualita- 
tively distinguish CNTs from the planar structures inves- 
tigated in [3, [13] ■ One can expect an essential role of 
these resonances in the formation of the CNT's thermal 
radiation. The black-body radiation distribution from 
multi-walled CNT bundles reported in Ref. ^] has been 
observed above the frequency range of surface plasmon 
modes and, what is more, is infiuenced by the nonhomo- 
geneous broadening due to the CNT length and radius 
dispersion and multi- walled effects. 

In this Letter we report calculations of the CNT ther- 
mal radiation spectra in the terahertz frequency range 
and demonstrate spatial coherence of the radiation. The 
model of the CNT optical properties based on the ef- 
fective boundary conditions jlSl ] and the fluctuation- 
dissipative theorem in the Callen-Welton form [3] have 
been utilized in our analysis. We show that the CNT 
thermal radiation is substantially different from the 
black-body radiation [3] and the radiation of planar 
structures with surface plasmons [3, [l3] ■ 

Consider an isolated single-wall CNT of the cross- 
sectional radius Rcn and the length L, aligned along the 
z axis of a cylindrical coordinate basis (p, </), z) with the 
origin in the CNT geometrical center. We restrict con- 
sideration to the case Rcn << 27r/fc, which implies the 
incident field to be slowly varied within the CNT cross- 
section; here k — uj/c, u is the electromagnetic field fre- 
quency, and c is the vacuum speed of light. The CNT is 
assumed to be placed into an optically transparent rar- 
efied medium; the CNT is in the thermal equilibrium 
with this medium. Optically, this medium is equivalent 
to vacuum. The Hamiltonian gauge, which implies the 
scalar potential to be equal to zero, is used for electro- 
magnetic field. Then, to calculate thermal fluctuations 
of the field in the presence of CNT we make use the 
fluctuation-dissipative theorem, which expresses the cor- 
relations of the vector potential in terms of the retarded 
Green tensor G(ri, r2, w) of a system (9|: 
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Im[G„m(ri,r2,w)] , 
(1) 

where n,m = p, (j), z, 0{lu, T) = huj/ [ex^(hjaj / ksT) — 1], 
T is the temperature, h and /cg are the Planck and Boltz- 
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mann constants, respectively. 

The equation G^°^(ri,r2, = (I + k-^Vr, 
VrJ G^°''(ri, r2, w) with I as the unit tensor and 
Vri as the operator dyadic acting on variables ri, re- 
lates the retarded free-space Green tensor with the 
corresponding scalar Green function G'*^°^(ri, r2, Ci^) = 
exp(iA;|ri — r2|)/|ri — r2|. For a given m, the column 
Glm(ri, r2, w) can formally be considered as a field vector 
induced in the point ri by a source located in the point 
r2. The scattering problem for that incident field is for- 
mulated by the equation (V^ x x — fc^) G(ri, r2, w) = 
47rl(5(ri — r2) imposed to the effective boundary condi- 
tions on the CNT surface We present the solution 
in the form of the simple-layer potential [i^l : 

G'„m(ri,r2, w) = Gl^^{ri,r2,uj) + '^^^^^ 

X j j(")(z;r2) j GW(ri,R,c^)d0dz, (2) 

-L/2 

where j^\z;Y2) is the normalized density of the axial 
current induced on the CNT surface by the incident elec- 
tric field Gim(R, r2,w), R = {Rcn,4>,z}- While deriv- 
ing ([2]) we assumed the incident field source distant from 
the CNT farther than its radius; therefore we can neglect 
the current ji™'' dependence on the azimuthal variable (f>. 

Eq. ([2]) with arbitrary ji'"'' satisfies the afore stated 
equation for the retarded Green tensor and the radiation 
condition at |ri — T2\ oo. Using the effective bound- 
ary conditions we obtain for jl™'(z;r2) the integral 
equation as follows 

rL/2 

/ ji"^^'; r2)/C(z - z')dz' + Cie-'^'- + Gae*'^'^ 

J-L/2 

^T- r G<f^{R',r2,io)dcj,'dz', (3) 

J-L/2 Jo 

where 

IC{z)^ (z) - r^—d^, (4) 

r = ^z^+ARl;j^ sin2(,/,/2), /Ci(z) = exp{ik\z\)/2ik, 
Gzzi^) is the CNT axial conductivity [l^, Gi and 
G2 are constants determined by the edge conditions 
ji™-' (±L/2; r2) — 5i|. The index m and the variable 
T2 appear in Eqs. ^ and ^ only as parameters. Note 
that these equations, as they couple the Green tensor of 
the system considered and the free-space Green tensor, 
play the role of the Dyson equation for CNTs. 

It is of importance that the role of scattering by CNT 
is not reduced to a small correction to the free-space 
Green tensor. This means that the iteration procedure 
conventionally used for the Dyson equation solving gets 



inapplicable in our case. Because of that, the direct nu- 
merical integration of Eq. ([3]) has been performed with 
integral operators approximated by quadrature formula 
and subsequent transition to a matrix equation. 

First term in square brackets in Eq. ([T]) is due to zero 
point energy of electromagnetic field and is further omit- 
ted. The second term comprises the CNT's thermal ra- 
diation and the black-body radiation of the surrounding 
medium. The latter one is represented by the correla- 
tor {A\^\yi)A^^'^*{y2))^ = D^nn}iri,r2,u;) of the black- 
body radiation vector potential A^-^K To separate off 
the CNT's intrinsic thermal radiation, we make use the 
method presented in Ref. Q (see the problems after Sect. 
77). The electric field intensity of the CNT's thermal ra- 
diation /aj(ro) — |E(ro)P is given by 

/.(ro) = fc2^[(K(ro)P)^-i?if,)(ro,ro,c.)]. (5) 

n=l 

Here, the relation En = —ikAn has been accounted for. 
Note that Eq. ([5]) also determines the intensity of the 
electric field fluctuations in the case when the CNT tem- 
perature is much higher than the temperature of the sur- 
rounding medium. 

By analogy with Eq. the vector potential An (ri) 
is written as 

L/2 27T 

+ ^ J j{z)jG^^)Ari,R,^)d<Pdz, (6) 

-L/2 

where j4^°''(r) is the vector potential of the free-space 
black-body radiation (i.e., of the black-body radiation in 
the absence of CNT). Second term in ^ describes scat- 
tering of the free-space black-body radiation by CNT. 
The density of the current j{z) induced on the CNT sur- 
face is found as a solution of Eq. ^ with iujA'^^^ (R, uj)/c 
substituted instead of Gim(R, r2, c^)). To calculate fim 
we utilize ^ and take into account that the correlator 
(AL'''(ri)A^'*(r2))^ is defined by Eq. ^ with the free- 
space Green tensor Gnm (ri , r2 , w) in the right-hand part. 

The thermal radiation spectra from (15,0) CNT at dif- 
ferent distances from the CNT axis are depicted in Fig. 
H^; Fig. [TJd presents one of the spectra in the logarithmic 
scale. For comparison, the black-body radiation inten- 
sity /i,^^(ro) = Auj'^e{uj,T)/c^ is also depicted. 

The spectrum depicted in Fig. [T^ demonstrates a num- 
ber of equidistant discrete spectral lines with decreasing 
intensities superimposed by the continuous background. 
Such a structure is inherent to spectra both in the far- 
field (dashed line) and near-field (solid line) zones. The 
peculiarity of the near-field zone is the presence of addi- 
tional spectral lines absent in the far-field zone. Thus, 
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FIG. 1: (a) Thermal radiation spectra of metallic (15,0) CNT 
in the cross-section zq = and at the distance from the CNT 
axis po — lOOL, (dashed line, left ordinate axis) and po = 0.5L 
(solid line, right ordinate axis). On the inset, the polariz- 
ability of the CNT is depicted, (b) Thermal radiation from 
CNT (solid line) and black-body radiation (dashed line) in 
the near-field zone. In both pictures L = 1 pm and T — 300 
K. The relaxation time r = 3 x 10~^^ s has been used under 
calculation of the axial conductivity a^z- 

the thermal radiation spectra presented in the figure are 
quaUtatively differ from both black-body radiation [§] 
and radiation of semi-infinite SiC samples [l^. In the 
latter case, the discrete spectrum is observed only in the 
near-field zone 

Comparison of the thermal radiation and the CNT's 
polarizability spectra depicted in Fig. [T^ reveals the coin- 
cidence in the far-field zone of the thermal radiation res- 
onances and the polarizability resonances. The latest are 
the dipole geometrical resonances of surface plasmons 0] 
defined by the condition Re[K(ijj)]L = 7r(2s — 1) with k{(jj) 
as the plasmon wavenumber; s is a positive integer. It 
should be noted that the polarizability (and the thermal 
radiation) resonances is found to be significantly shifted 
to the red as compared to the perfectly conducting wire 
of the same length because of the strong slowing^own 
of surface plasmons in CNTs: Re[K(a;)]/fc « 100 [l8[. In 
particular, for L = I /im, the geometrical resonances fall 
into the terahertz frequency range. The attenuation is 
small in a wide frequency range below the interband tran- 
sitions. Additional spectral lines in the near-field zone 
are described by the condition Re[K(ti;)]L = 27rs and are 
due to the quadrupole geometric resonances. Thus, the 
resonant structure of the thermal radiation spectra is en- 
tirely determined by the finite-length effects. Note that 
similar structure of the thermal radiation spectra is pre- 



dicted for 2D-electron gas [24|. Resonances in Ref. (221 
are due to excitations of other physical nature — optical- 
phonon modes of the barrier material. 

The presence of singled out resonances illustrated by 
Fig. [T^ allows us to propose metallic CNTs as far-field 
and near-field thermal antennas of a new type (thermal 
antennas based on photonic crystals have recently been 
considered in Ref. 0). Taking into account the high 
temperature stability of CNTs, the CNT thermal anten- 
nas can be excited by Joule heating from the direct (low- 
frequency) electric current. 

Accordingly to Refs. [1, 0, Q the maximal efficiency of 
vibrator CNT-antennas is reached at frequencies of the 
surface-plasmon dipole resonances. Figure [1^ shows that 
the intensities of spectral lines of the thermal radiation go 
down with the resonance number much slowly than the 
polarizability peaks. This means that the signal/noise 
ratio for the CNT-based antennas is maximal for the first 
resonance and decreases fast with the resonance number. 
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FIG. 2: Same as Fig. [T] but for the (23,0) semiconducting 
zigzag CNT. 

As different from metallic CNTs, semiconducting ones 
do not reveal isolated resonances in both far-field and 
near-field zones, see Fig. [5] as an illustration. Such a pe- 
culiarity can easily be understood by accounting for the 
strong attenuation of surface plasmons in semiconducting 
CNTs whereas the slowing-down remains of the same or- 
der. That is why in this case the Q-factor of geometrical 
resonances turns out to be substantially smaller and they 
do not manifest themselves as separated spectral lines. 
In the same way, the thermal radiation intensity of semi- 
conducting CNTs is substantially smaller than that of 
metallic ones and displays qualitatively different spectral 
properties in the near-field zone: monotonous growth of 
the intensity with frequency inherent to metallic CNTs 
changes into monotonous declining (see Fig. [2^). As the 
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thermal spectra are strongly dependent on the CNT con- 
ductivity type and length, the near-field thermal radia- 
tion spectroscopy proposed in [l^ for testing the surface- 
plasmon structures, can be expanded to CNTs. 

Figure [^b demonstrates that in the frequency range 
considered the black-body radiation intensity consider- 
ably exceeds the thermal radiation of semiconducting 
CNTs: <C /i . In the vicinity of geometrical res- 
onances the same property is inherent to metallic CNTs, 
see Fig. [Ud. This means that CNTs as building blocks 
for nanoelectronics and nanosensorics possess uniquely 
low thermal noise and, thus, provide high electromag- 
netic compatibility on the nanoscale: Their contribution 
to the electromagnetic fluctuations in nanocircuits is neg- 
ligibly small as compared to the contribution of dielectric 
substrate. More generally, the later example illustrates 
peculiarity of the electromagnetic compatibility problem 
on the nanoscale motivating future research investments 
into the problem. 

In vie of the CNTs' application as tips for near-field op- 
tical microscopy , we have studied the spatial structure 
of the electromagnetic fluctuations near CNT, character- 
ized by the normalized flrst-order correlation tensor 

(A„(ri)A;,(r2)). 



3m«(l"l>I"2,t^) 



V(|A.(ri)P)^(|An(r2)P). 



(7) 



Calculations of the axial-axial component of this tensor 
in the near-fleld zone are presented in Fig. [31 Because of 
the strong slowing-down of surface plasmons in CNTs, 
the vacuum wavelength of the dipole geometrical reso- 
nance A ~ lOOL exceeds signiflcantly the corresponding 
coherence length, which is estimated from the peak width 
in Fig. [3] (solid line) to be about lOL. At the quadrupole 
resonance frequency (dotted line) and far away from res- 
onances (dashed line) the coherence length considerably 
increases. Physical interpretation of such a behavior is 
analogous to that given in Ref. [3] for planar struc- 
tures: The coherence length is of the order of magnitude 
of the spatial scale of the dominant mode. At the dipole 
resonance the dominant mode is surface plasmon those 
spatial scale is dictated by ^ 2'kX. Far away from 
the resonance the contribution of volume modes consid- 
erably increases extending the spatial scale to ~ A/2. 

Concluding, the thermal radiation from single-wall 
flnite-length CNTs in the terahertz range has been the- 
oretically investigated. Discrete spectral lines in the ra- 
diation spectra of metallic CNTs, originated from the 
geometrical antenna resonances, is predicted to exist in 
both near- and far-field zones. The effect allows a con- 
cept of metallic CNT as a thermal antenna and is of 
importance for the CNT-spectroscopy, the nanoantenna 
design, the high-resolution near-field optical microscopy 
and the thermal noise control in nanocircuits. 
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FIG. 3: First-order correlation function g'^^^} {v\,V2,oj) for 
metallic (15,0) CNT. n = {0.5L, r2 = {0.5L, 0, 0}, and 

(j) has an arbitrary value. 
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